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Abstract 

The capillary instability of liquid crystalline (LC) jets is considered in the framework 
of linear hydrodynamics of uniaxial nematic LC. The free boundary conditions of the 
problem are formulated in terms of mean surface curvature H and Gaussian surface 
curvature Q. The static version of capillary instability is shown to depend on the elas- 
ticity modulus K, surface tension do, and radius tq of the LC jet, as expressed by the 
characteristic parameter x = K/a^rQ. The problem of capillary instability in LC jets is 
solved exactly and a dispersion relation, which reflects the effect of elasticity, is derived. 
It is shown that increase of the elasticity modulus results in a decrease of both the cut 
off wavenumber k and the disturbance growth rate s. This implies enhanced stability 
of LC jets, compared to ordinary liquids. In the specific case, where the hydrodynamic 
and orientational LC modes can be decoupled, the dispersion equation is given in closed 
form. 



Key words: Jet, Rayleigh Instability, Nematic Liquid Crystal 



e-mail: lfel@techunix.technion.ac.il 



1 Introduction 



The breakup of liquid jets, that are injected through a circular nozzle into stagnant fluids, 
has been the subject of widespread research over the years. Previous studies that followed 
the seminal works of Lord Rayleigh have established that the complex jet flow is influenced 
by a large number of parameters. These include nozzle internal flow effects, the jet velocity 
profile and the physical state of both liquid and gas. Notwithstanding the fact that the 
hydrodynamic equations are nonlinear, the linear stability theory can provide qualitative 
descriptions of breakup phenomena and predict the existence of different breakup regimes. 

Rayleigh showed [1] by using a linear theory that the jet breakup is a consequence of 
hydrodynamic instability, or more exactly capillary instability. Neglecting the effect of the 
ambient fluid, the viscosity of the jet liquid, and gravity, he demonstrated that a cylindrical 
liquid jet is unstable with respect to disturbances characterized by wavelengths larger than 
the jet circumference. He also considered the case of a viscous jet in an inviscid gas and an 
inviscid gas jet in an inviscid liquid [2]. Weber [3] generalized Rayleigh's result for the case of a 
Newtonian viscous liquid and showed that the viscosity tends to reduce the breakage rate and 
increase the drop size. Chandrasekhchar [4] considered the effect of a uniform magnetic field on 
the capillary instability of a liquid jet. A mechanism of bending disturbances and of buckling, 
slowly moving, highly viscous jets, was presented by Taylor [5]. Further developments of the 
theory in Newtonian liquids was concerned with additional factors such as the dynamic action 
of the ambient gas (leading to atomization of the jet), the nonlinear interaction of growing 
modes that lead to satellite drop formation, and the spatial character of instability (see [6], 
[7])- 

The capillary instability in jets, comprised of non-Newtonian suspensions and emulsions, 
presents a different category of cases which are governed by power-law (pseudoplastic and 
dilatant) liquids. The effective viscosity of the pseudoplastic liquid decreases with growth of 
strain rate, whereas in dilatant liquids, it increases [7]. The behaviour of capillary jets of 
dilute and concentrated polymer solutions suggest a strong influence of the macromolecular 
coils on their flow patterns [7]. Free jets of polymeric liquids, that exhibit oscillations, are 
reported in [8]. 

Recently the idea of Rayleigh instability was applied to tubular membranes in dilute ly- 



otropic phases [9]. Their relaxation, following optical excitation, is characterized by a long 
time, and can be described by means of hydrodynamic approach [10]. Bending deformations 
of such membranes are governed by the Helfrich energy [11] which depends on the curvature 
of the tube. Thus, competition between the surface tension and curvature energy of the wa- 
ter immersed membrane renders the initial shape of the tube unstable. The hydrodynamic 
formalism used in [10] and the hydrodynamics of fluids with inner order such as liquid crystal 
(LC) [12] have similar features. In [10] the order parameter stands for a unit vector normal to 
the membrane surface. In contrast, the order parameter Q of a LC fluid, is defined throughout 
the space it occupies. 

The continuum theory of LC phases has emerged as a rigorous part of condensed matter 
theory. The hydrodynamics of the LC phases was developed during the 70-80th and its 
predictions were successfully confirmed in many experimental observations. The combination 
of viscous and elastic properties is likely to produce new evolution patterns of hydrodynamic 
instabilities, in the context of Benard-Rayleigh, Marangoni and electrohydrodynamic effects 
[13], which cannot occur in ordinary liquids. However, its capillary instability, when in the 
form of a jet, was not considered as yet. In particular we refer to the uniaxial Hematic phase. 

The instability of a LC jet poses an additional challenge with respect to the effects listed 
above. This applies already within the framework of linear stability theory. The LC class of 
fluids seems to provide a good example of unique properties, as compared to polymer solutions. 
The elastic properties of a LC are expected to change the evolution patterns of jets which are 
made from them. In this work we derive a rigorous mathematical model of capillary instability 
for isothermal incompressible nematic LC jets. This model shows how the combined viscous 
and elastic properties of LC fluids determine the boundary conditions at the free surface, and 
the range where instability prevails. 

2 Hydrodynamics of a liquid crystalline jet 

In this Section, we formulate first the problem of capillary instability and then derive the 
basic equations which govern the linear hydrodynamics of a liquid crystalline jet. The flow 
of a nematic LC is described by a set of differential equations supplemented by boundary 
conditions on the LC free surface: continuity equation, Navier-Stokes equation of visco-elastic 



LC, and Lesli-Ericksen equation of angular motion of the director n(r, t). 

The basic notations and linear hydrodynamic equations of uniaxial nematic liquid crystals 
follow the theory given in [12], [14], [15]. 

2.1 Basic notations and variables 

The following basic variables describe the nematic LC medium: velocity V(r, t), pressure 
P(r,t) and LC-director n(r, t). The initial values of the functions will be denoted by "o", 
either as a subscript or superscript. The following notations, which are commonly accepted 
in the theory of LCs, are used henceforth: 

1. The free energy density E d of deformed non-chiral uniaxial nematic LC, given in 
quadratic approximation in terms of the derivatives dn/dxj reads 

2E d = K x div 2 n + K 2 (n, rot n) 2 + K 3 [n x rot n] 2 , (1) 

where (a, b) and [a x b] denote scalar and vector products of vectors, and Ki > 0, % — 1, 2, 3 
are known as the Frank elasticity moduli. In the vicinity of a phase transition oc Q 2 [13] 
and in the isotropic phase they vanish. 

2. The bulk molecular field F and the Ericksen elastic stress tensor r ki) which set the 
equilibrium distribution of the n-field in a LC, are determined by the following variational 
derivatives 1 

F = M - n(n, M) , or F i = (<% - mn^Mj , 

where 

A 9Ed 9Ed .= 9Ed f) A 

dx k d{d k ni) drii o{d k ni) dx k 

i.e. 

M = Ki grad div n — K 2 {(n, rot n)rot n + rot ((n, rot n)n)} + 

i^ 3 {rot [n x [n x rot n]] + [[n x rot n] x rot n]} , 

T ki = Ki fediv n + K 2 (n, rot n)n m e mki + K 3 [[n x rot n] x n] m e mki , (4) 
1 Here and throughout the paper, unless noted otherwise, we apply the summation rule over indices which 
are repeated in a tensor product, e.g. atjbjh — J2j a ijbjk- 



(2) 
(3) 



e m ki is a completely antisymmetric unit tensor of the 3rd rank (Levi-Civita tensor). 

3. If the deviations of the director n = n° + n 1 from its initial orientation n° are small, 
then 



n o _..u 



n° y = 0, n° z = 1 , 1 » nl,n\ » n\ ~ (n x f , (nj) 2 , (5) 



and simple algebra yields the following linear approximation 

- d 2 n l ~ d 2 n x 8 

M x = Kn 1 x + (K 1 -K 2 )^, M y = Knl + (K 1 -K 2 )^, M z = (K, - K 3 ) -div n 1 , (6) 

where K = + K 2 j^ + x 3 |L and by virtue of (2), F x = M x , F y = M y , F z = 0. 

If further simplification through single elastic approximation K\ = K 2 = K 3 = K is applied, 
then 

r^2 r^2 ^2 

F x = KA 3 nl, F y = KA 3 n y , F z = , A 3 = ^ + ^ + ^ , (7) 

where A3 is the three-dimensional Laplacian. Similar considerations regarding the Ericksen 
stress tensor r ki give 



T~xx — T yy ~ r zz ~ K\ div n 1 , T xy — Ty X — K 2 



dx dy 



The stresses given by (8) do not contribute to the non-dissipative stress tensor T- k used in 
the linear hydrodynamics of LCs (see (9) below). 

4. The reactive (non-dissipative) T\ k and dissipative Tf k stress tensors are defined as follows 

On • A i 
T\ k = -P 5 lk - r kj - - ( ni F k + n k Fi) + - (n t F k - n k Fi) , (9) 

Tf k = 2rnT ik + (772 - 771) 5 ifc divV + (771 - 772 + r]i)(SiknjT jm n m + mn k <liv\) + (10) 

(773 - 2771) (niT kj nj + rikTijUj) + (771 + 772 + 775 - 2773 - 2ri A )n i n k n j n m T jm , 

where the antisymmetric VL ik (vorticity) and symmetric T ik parts of the derivative d k Vi read 

_lfdV k OVA _lfdV k dVA 



2 y (9xj dx k ) ' 2 y <9x j ctefc , 

Five independent viscous moduli 77^ , kinetic coefficient A, and rotational viscosity 71, deter- 
mine the dissipative stress tensor Tf k , the 4th-rank viscosity tensor r) ik j m , and the dissipative 



function D in the absence of heat fluxes 

D Vikjm\ ik^- ' jm F , T^ k Vikjm^- jm i (^^) 

ll 

Vikjm = Vl{iijikm + CkjCim) + i.V2 ~ Vl)£,ik£,jm + ^(?V^£fcm + ^k^jiim + ^i^m.£fcj + 

+n fc n m ^) + i] A (nin k ^ jm + njn m £ ik ) + V5nin k njn m , £ ik = 5 ik - rnn k . 

The tensor rj ik j m consists of five independent uniaxial invariants [12] and is highly symmetrical 
Vikjm = Vkimj — Vjmik- The requirement that D be positive translates into, 

77i > , ?72 > , r/ 3 >0, 775 > , mV5>vl, 7i >0. (13) 

The parameter A is close to +1 or —1 for rod-like or disk-like molecules, respectively. If the 
liquid is visco-isotropic, then A = 0. 

5. The hydrodynamic reactive (non-dissipative) m r and dissipative m d fields are defined 
as follow 

m'i = -(V, V 3 )rii + n k Q ki + A ^jT jk n k , m d = — F , (14) 

7i 

where V3 is the three-dimensional gradient operator, (V3) 2 = A 3 . 

6. The surface tension a of nematic LC is given by [16], 

a = (7 + 0-1 (n, e) 2 , (15) 

where o"o and cri are isotropic and anisotropic surface tension moduli respectively, and e is a 
unit normal vector to the LC surface. 

7. In the case of an incompressible LC (r/'" = r/'^, vT = 0) the tensors Tf k and Vikjm take 
the following form, 

T* n = 2<T lfc + (r? 3 n - 2<) (n,^ (16) 

V3 

Vikjm = Vl (^ij^km + ^kj^irn) H ^~( n i n j^km + n k n j£im + n i n m^kj + n k n m£ij) + % ^i^fc^ j n m j 

where "in" denotes incompressibility condition. 

8. The first two terms in (10) and the second equation of (12) correspond to ordinary 
compressible liquids with isotropic invariance. The simplification of (10) results from 



Vs = 2?7i , v\ = v\ - v[ , Vs = v\ + v[ , 



so that 



Tik = 277iT ifc + (r? 2 L - r][) 5 ik divV , rfc kjm = r}[($ijhm + S im 5 kj ) + (% L - r][)5 ik 5 jm . 

The coefficients r][ and rfe — r^J- are known as the first and second isotropic viscosities. 

9. Another system of viscous moduli ai (called Lesli viscosities) relate dissipative and 
kinetic moduli in the following way 2 

«4 x 72 

T]i — — , A = , f] 5 = ai + 0(4 + a 5 + a 6 , 71 = a 3 - a 2 , 72 = «3 + « 2 , 



with the support of Onzager-Parodi relation [17] 0:3 + «2 = «6 — «5- In the vicinity of phase 
transition, the viscous moduli ctj have different dependences upon the order parameter Q: 
CKi oc Q 2 , CK2, «3, 0:5, a6 oc Q, 0:4 oc Q° [13]. 

Tables 1 and 2 (see Appendix) summarize viscosities and other physical parameters that 
characterize the most frequently used and well studied nematic LC, also known as MBBA and 
PAA. 

2.2 Basic equations 

The complete system of hydrodynamic equations for nematic LC reflect the conservation laws 
of mass, and of linear and angular momenta. 
1. Continuity equation 




(17) 



dp 

dt 



+ div (pV) = . 



(18) 



2. Navier-Stokes equation for visco-elastic LC 




(19) 



dx k 



3. Lesli-Ericksen equation of angular motion of the director n(r, t) 




(20) 



2 The correct expression for r] 5 is given in [14]. 



The last equation is written for a negligible specific angular moment of inertia Jlc of the LC, 
namely, Jlc <C p r$, where r is a characteristic size of the system. This is true in our case, 
where r denotes radius of the jet. 

Consider an isothermal incompressible jet, flowing along the z axis, out of a nozzle at a 
velocity V. The initial orientation of director n° is assumed collinear with V. The deviation 
from initial values of the director and pressure are defined as n 1 = n — n°, and Pi = P — P 
respectively, where P = cr/r is the unperturbed pressure within the cylindrical jet. Applying 
the linear approximation In 1 1 <C 1, equations (18)-(20) are simplified as follows 



dVi 

divV = 0, 

dn} 



dPi 

dxi ' dxk 



+ 



dT«? 



+ 



— n° divF — 



(n°,V 3 )F,, 



n° k ttki + A ^ijTjknl + —Fi , f 9. = ^ - nM , i, j, k = x,y,z 



dt - - - 7i 

Choosing n° z = 1, gives F z = and hence 



(21) 



= 



dV r . dV„ dV x 



+ 



dV x 



p- 



dx dy 
dPi 



+ 



dz 



dt 

dVy 

dV z 
dt 

<K 

dt 



dx 

dPi 

dy 
_dPi 

dz 
A + 1 dV x 

2 dz 
A + 1 dV y 



+ 
+ 
+ 



Pi A 2 + /3 2 
fa A 2 + p 3 



d^_ 

dz 2 

dz 2 
d 2 



V x + {(3 2 -Pi) 

Vy + (& - f3l) 



d 2 V z _ A + 1 dF x 
dxdz 

d 2 v z 



dydz 



2 dz 
A + 1 dF y 
2 dz 



(22) 
(23) 



dz 2 



V, 



A 



divF , 



+ 



A - 1 dV z 



dx ^ 71 



(24) 



+ 



A - 1 dV, 



+ 



dn\ 



0, 



2 dz 2 dy 7l ' 9t 

where A 2 = J^- + Jp- is the two-dimensional Laplacian, = r^ n , /3 2 = ^3 n /2, /?3 = i]'" 

and F x ,F y are given in (7). As isotropic viscosity means $ = /3, the above mentioned liquid 

crystals, MBBA and PAA, are clearly far from being isotropic (see Tables 1, 2 in Appendix). 

In order to make the problem more specific and easier to solve, we consider axisymmetrical 
disturbances in a system of cylindrical LC jet, with radius r and subject to the single elastic 
approximation [Ki = K). This provides the simplest approximation which still preserves the 
influence of elastic forces, on the hydrodynamics of an incompressible and elastic LC. In this 



case 



P 
P 

7r 





dV r 

dt 

dt 
dr^_ 
dt 



dV z dV r V r 
— t _| + — 

dz dr r 



dP 1 
dr 
dP 1 



dz 



Pi A 



*2c 



dz 2 



v r + (p 2 -Pi) 



d 2 v z 

drdz 



Pi 



dF r 

dz 



dz 2 



liPr 



dV r 
dz 



/3 2 A 2c + 



V z - p 2 



dF r 



dr 



El 

r 



0, 



where 
A 2c 



d 2 1 d 
1 F 



K { A 2c - 4 + 



dz 2 ] Ur ' 



Pi 



A + l 



P2 



A - 1 



(25) 
(26) 

(27) 

(28) 

(29) 



Uz- ! 2 2 

Equations (25)-(28) describe ordinary linear hydrodynamic behaviour of isotropic incom- 
pressible liquids provided that the LC properties vanish: K, ^ — > and fa = f3. The result is 
the well known continuity and linearized Navier-Stockes equations. 

<9V 



divV = , p- 



dt 



-VPi + /9A 3 V . 



(30) 



2.3 Boundary conditions at free surface 

Boundary conditions at the free surface of a liquid crystal state that the jump in normal 
stress consists of two parts: one depends on the surface tension a, and the other on the elastic 
disturbance W e i ast of the uniform director field n (r). Assuming that no tangential stresses 
exist at the free surface, the boundary conditions can be expressed as, 



do" 

(T^ k + Ti n )e k + (2aH + W elast )e i + -^ ^0 at r = r , 



(31) 



where are the components of the normal unit vector e in the reference frame of the LC- 
cylinder, and H, = 1/2 (1/Ri + I/R2) denotes mean surface curvature with principal radii Ri 
and R 2 . 

The non-hydrodynamic part of boundary conditions at the free surface holds, provided 

that the scale of deformation of the initial surface is considerably larger compared to the 

molecular length of LCs 3 . This dictates a tangential behaviour of a smoothly disturbed 
3 Strictly speaking, this assumption is correct when an equilibrium distribution of director field n(r) is free 
of singularities. The problem of minimal surface of LC drop presents another situation wherein an essential 
rearrangement of the field n(r), at the surface, can diminish the total energy by destroying the disclination 
core within the drop. 



director n at the free surface, e 2 < e r ~ 1 : 

(e, n) = -> e z + nl = at r = r . (32) 

The last constraint cancels the gradient term in (31). Finally we come to the boundary 
conditions in the linear approximation of the variables nl,V r ,V z , Pi 

T; r + T r d l" + 2aH + W elast = 0, T; r + T z d r = 0. (33) 

Substitution of the expressions for the reactive and dissipative stress tensors gives 

2p 1 T rr - Pi= 2a (Ho — H) — W elast , 2[3 2 T zr = fi 2 F r at r = r . (34) 

where Ho = l/2r is the initial mean curvature of the LC-cylinder. The equations for a jet 
surface, disturbed by a wave ((z,t), and its radial velocity d(/dt, are given by 

dC 

r(z, t) =r + ((z, t) , V r = — at r = r , (35) 

where ( r is the radial displacement of a surface point. The principal radii of the surface 
curvature, in the context of linear approximation with respect to (, and its derivatives can be 
expressed as, 

1 1 1 C 1 d 2 ( 

This transforms the boundary conditions (32), (34) into 





dz 


V r 


d( 
" dt ' 




(37) 


2(3 2 T zr 




i 






(38) 


Pi — 2(3iT rr 


= -o"o ( 


'C 


dz 2 J 


+ w elast . 


(39) 



The term W e i as t deserves further discussion. It reflects the existence of normal stresses, at the 
surface, which arise due to the resistance of the uniformly orientated continuos LC media to 
a surface disturbance. W e \ ast vanishes in undisturbed LC jets and it depends linearly on the 
elastic modulus K, radius r and derivatives of (. Moreover, an invariance of the problem 
with respect to inversion of the z-axis requires sole dependence on derivatives of even orders. 
An explicit expression for W e i ast is derived in Section 3.1. 



3 Plateau instability in a LC cylinder 

Before proceeding to tackle the sophisticated mathematics of equations (25)-(28), as supple- 
mented by boundary conditions (37)-(39), capillary instability of the LC cylinder is discussed. 
This is done by applying the Plateau considerations [18] on the figures of a liquid mass with- 
drawn from the action of gravity. 

Consider a LC cylinder with a surface disturbed as specified by (35), where ( = Co c °s kz, 
Co is small compared to r , and k = 2n/A, A being the disturbance wavelength. The idea of 
Plateau, applied here, is to find such cut-off wavelength A s of the disturbance, that defines 
breakage of the cylinder into droplets with due decrease of the total energy. 

The volume v enclosed within one wave length is given by 



(40) 



where r in the right h.s. of (40) is given as a second order expansion ( . The total energy E 
of the LC cylinder with a disturbed director field n(r) is given by 

E = (T [ ds + ^- f (div 2 n + rot 2 n) alv . (41) 
J s 2 J ' v 

The static director field n(r) can be found from equation (29) and the attendant boundary 
condition (37) 

/ 1 d 2 \ dC 

n°=l, F r = -> (A 2c -- + — Jnj = 0, n\ = A at r = r . (42) 

Equation (42) is satisfied by the following solution, which is finite at r = 

nl(r,z) = - k< *° . h(kr) sin kz , (43) 
h{kr ) 

where I m {x) is a modified Bessel function of order m. The contribution of elastic forces is 
determined by 

Ko 



div 2 n + rot 2 n = k 2 
where 



h(kr ) 



[A\{kr) sin 2 kz + A\(kr) cos 2 kz] (44) 



My) = ^- + -h(y), A 2 (y) = h(y). 
dy y 



A simple integration of (41) gives 

£ = 2wo U + \k 2 Q\+-K 



\ k( 


7 


[h(kr )_ 


Jo 



i>kr 

[A\{y) + Al{y)]ydy. (45) 



Inserting tq from (40) into the first term above, we obtain 



-Co 2 



S - 2a ^v~ = a ^ (w 2 -l) + ^K 
Zr Z 



C m 


7 




Jo 



[Al(y)+A*(y)]ydy, w = kr . (46) 



The positive root w s = k s rn, given in the right h.s. of (46), determines the cut-off wavelength 
A s of capillary disturbances, which renders the LC cylinder unstable. Subsequent disintegra- 
tion into detached masses is favored by the decrease in S 



K - 1) + *75 
1 i 



[Aj(y) + A 2 2 (y)} ydy = , 



K 



x 



(47) 



\vj s ) .In a r 
where the subscript " s" denotes the static nature of Plateau instability. 

The quadratic approximation (1) with respect to the derivatives dn/dxj, which provides 
the basis for the Frank theory, makes the expression (47) correct only in terms of the w 2 s 
approximation. Indeed, the power of w s in (47) should not exceed 2, otherwise the calculation 
becomes inconsistent. Thus, we get 



ir( 2 1 
8 — 2cr V / v™ = a — - (w 2 — l) + TrKk 2 ^ and w s = ■. 

2r VI + 2x 

The asymptotic behaviour of vj s {>c) shows two important limits: 



(4* 



w. 



x if x<l: w s 



1 - — ] if x > 1 . 

4x . 



(49) 



Figure 1 shows a plot of k s rn vs. x for Plateau instabilities in LC and in ordinary liquid. 
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Figure 1: Universal plots of k s r vs. x for Plateau instabilities in LC cylinder (plain line), 
and in ordinary liquid k s r = 1 (dashed line). 



The corresponding asymptotic cut-off wavelength A s are obtained as 



A s = 27rr (1 + x) if x<l; A, = 27r \/^ v^o U + ^ ) if x>l. (50) 



This result shows that k > k s increases the total energy £ of the disturbed system, whereas 
k < k s decreases it. According to (49), there are two marginal regimes of instability: 

• Capillary regime ro ^> Kjo§. Here A s is close in value to the circumference of the 



This regime must apply to a wide range of nematic LC, since the common values of K ~ 
1CT 11 J/m [13] and a ~ 1(T 2 J/m 2 [20] lead to K/a ~ 10~ 9 m. This value is evidently 
smaller than the presently attainable radii of the jet. 

• Elastic regime r <C K/<j . This case reflects the dominance of elastic deformation and 
predicts an unusual behaviour for A s ~ y/r~o. 

This regime cannot be reached by simple increase of the elastic moduli since their magnitude 
is determined by K ~ nT/a, where kT ~ 10~ 20 J is the Bolzmann thermal energy at room 
temperature, and a ~ 10~ 9 m denotes molecular length of LC. In contrast, the effect of surface 
tension can be diminished by surfractants or by charging the surface of the liquid. In the latter 
case the charge can virtually eliminate the effect of surface tension and provide the conditions 
where the elastic forces predominate. 

3.1 W e i as t and Gaussian surface curvature 

The straightforward way to derive an expression for W e i ast is to solve the elastic problem for 
the stresses existing on a deformed axisymmetric surface of a LC cylinder. This relates to the 
Plateau instability, which obviates the need to repeat the entire procedure. 

When we turn from Plateau considerations on the static instability of LC cylinders to the 
capillary instability of LC jets, the question is whether the cut-off wavelengths of both the 
static A s and hydrodynamic problems coincide. This question was skipped by Rayleigh 
in his studies of isotropic viscous liquids, since for ordinary liquids both cut-off wavelengths 
always coincide A s = A d . This identity reflects a deep equivalence principle of the bifurcation 
point for non-trivial steady state of dynamic system, and the threshold of static instability 
concerned with a minimum of its free energy £ [4]. 

Making use of A s = A d we construct the term W e i ast which enters the boundary condition 
(39). To this end, we examine and represent the total energy (48) as follows 



cylinder and the elastic deformation contribution f E d dv, to the total energy £, is negligible. 




(51) 



Next, we compare the expression within the brackets with the right h.s. of (39). This gives 
W e i ast , which generates the elastic contribution in (51) 

W^ = 2K 5 , C = ^ = -Ig, (52) 

where Q is the Gaussian surface curvature in accordance with (36). Thus the final expression 
for boundary conditions (31) is based on two fundamental invariants of the surface curvature, 
i.e. mean surface curvature H, and Gaussian surface curvature Q. 



4 Dispersion relation 



Rayleigh was the first to observe [1] that contrary to Plateau, the instability problem is not so 
definite. The mode whereby a system deviates from unstable equilibrium must depend on the 
nature and characteristics of the small displacements to which this system is subjected. In the 
absence of such displacement, any system, however unstable, cannot depart from equilibrium. 
These characteristics, being hydrodynamic, reflect the effect of viscosity, which predominates 
over that of inertia. In the case of ordinary liquids, the mode of maximum instability which 
corresponds to the wavelength A R = 4.508 x 2r exceeds the circumference of the liquid 
cylinder. We anticipate that the instability of LC jets possesses similar features. 

The fact that a velocity potential does not exist in an anisotropic visco-elastic liquid, 
dictates a standard approach to this problem which was elaborated first by Rayleigh [2]. Let 
us define the Stokes stream function ^(r,t) and a director potential Q(r,t) as, 



r dz 



19* ! oe 

V z = -— and n r = — , 
r or or 



(53) 



so that the continuity equation (25) holds. From the other three equations (26)-(28) we have 



dP l 



ld^ 



dz 

o 2 e 



ld_ 

r dr 
1 



for 



drdz \r dr 

d (\m 



dr \r dr 



ld_ 

r dz 

d 2 m 



dz 2 



' r _d_ /lcNA | d 2 ^ 
dr \r dr J dz 2 



dV 



P ~dt + ^ 



dt 



d (\m 



^ 2 dr \r dr ) pl dz 2 



d 2 m 



drdt r 
Applying the commutation rules give, 



1 



+ -F r , 

7i 



F r = K A 2c + 



d^_ 

dz 2 



i_\ d& 

r 2 J dr 



(54) 
(55) 
(56) 



which facilities simplification of the above equations. Assuming that an axisymmetrical dis- 
turbance, characterized by a wavelength 2n/k, increases exponentially in time with the growth 
rate s, gives, 



st+ikz 



0, C, Pi, F r } = {i 4>(r), i 6{r), q(r), p(r), i f(r)} x e 



Inserting (57) into (54)-(56) gives rise to the following amplitude equations 



(57) 



1 dp 



d ( 1 dip 



1> 



kp 

86 
dr 



ld_ 

r dr 



or \r or I r 



Or \r or J r 71 or 



(58) 
(59) 
(60) 



The new variables in (57) require reformulation of the boundary conditions (37)-(39) as 
follows, 



1. 99 J 
kq = — , sq — k— 

or r 



H2 d (ldip\ 2 ip 



dr \r dr 



(t^ + ' T (61) 



where 



r = (7 U 



' l) +2 K-k 2 . 



The real form of the amplitude equations (58)-(60) and boundary conditions (61) imply that 
(57) divides the five variables into two groups: Pi,( and \&,0,F r . These groups are shifted 
with respect to each other by the phase angle n/2. 

4.1 Reduction of the amplitude equations 

In this Section we perform a standard procedure for the simplification of the amplitude equa- 
tions (58)-(60). Substituting / from (60), into the other amplitude equations we get 



I dp d (ld^\ 2 ip dd 

k or or \r or I r or 



/ 1 9 

kp = -— < r 

r or 



or \r or I r 



1 d ( dO 







,2^ K 9 



d fldip 
' dr I r dr 



71 dr 



Id ( dd 



r dr \ dr 



r or \ or 



(62) 
(63) 
(64) 



where 



B x = fa - 71//1//2 , B 2 = (3 2 + , B 3 = (3 2 + 7iA«2 > B 4 = (3 3 - 71//1//2 • 



(65) 



and B 2 > 0, B 3 > by virtue of (13). Let a new stream function x be defined as ip = r dx/dr. 
The orientational •& and kinematic v , viscosities, as well as other auxiliary functions, are defined 
by the following relations 

#=-, i* = ^ , u 2 = k 2 + -, w 2 = k 2 + - Q , - « 1 - u 2 < w 2 , (66) 
7i P Vi v Vi 

where the first inequality in (66) applies to known nematic LC fluids (see Tables 1, 2 in 
Appendix). Using the new notations we find the first integrals of the amplitude equations as, 

V - = (B 1 A 2c -B 2 u 2 2 ) X + s^ 1 e, (67) 
kp = (B 3 A 2c - B 4 ul) A 2cX - s7iA*2A 2c , (68) 
= (/i 2 A 2c + fiik 2 ) x + $ (A 2c - w 2 ) . (69) 

Next, we eliminate the pressure amplitude p from (67) and (68). This gives, 

[B 3 A 2 2c - (B^ 2 + B A ul)A 2c + B 2 u 2 2 k 2 ] X -s ll (fi 2 A 2c + fi 1 k 2 )e = 0, (70) 

{li 2 A 2c + ii 1 k 2 )x + '&{A 2c -w 2 )e = 0. (71) 

Diagonalizing a matrix of operators in (70) and (71) we obtain the following homogeneous 
equations for the functions x( r ) an d #(r), 

[D 3 A\ C - D 2 A\ C + D x A 2c _£) ]^j = ^j , ( 72 ) 

where 

D = k 2 (<dB 2 u 2 2 w 2 - s^fijk 2 ) ,£>i = # (B^w 2 + B 2 k 2 u\ + B 4 iu 2 u2) + 2s- fl fi lf i 2 k 2 , 
D 2 = t? (#iA; 2 + 5 3 w 2 + B^ul) - s 7l /i 2 , D 3 = $B 3 . (73) 

It is easy to verify that all coefficients Dj are positive, if the conditions that all Bi > 
and /i 2 <C 1, <C 1 are satisfied. The latter are in a good agreement with numerous 
observations in nematic LCs [13]. 

Further factorization (recalling that D 3 > 0) of the polynomial differential operator gives 

D 3 A\ C - D 2 A\ C + D,A 2c -D = D 3 (A 2c - m\) (A 2c - m 2 ) (A 2c - m 2 ) . (74) 

Equation (74) facilitates the following finite solutions of equations (72) 

^ C ^ G 

X(r) = Y] —Io( m 3 r ) > #( r ) = J] — I o( m j r ) > ( 75 ) 



where the second fundamental solutions that diverge at r = were excluded, Cj and G 3 are 
indeterminate coefficients and m? are three generic 4 roots of the following cubic equation 

3 3 3 

3=1 6 j+k 6 j=l 6 

The coefficients Gj can be expressed through Cj, once (75) is inserted into (71), 

1 U\k 2 + « 2 m^ 

G^^gjCj, 93= w2 _ m] \ J = 1,2,3. (77) 

The amplitude of the pressure p(r), the stream function ip(r) and the displacement of a point 
on the surface <r(r ) are easily found from (60), (67), (71) and (77), 

3 I- s 
p(r) = — C j I o( m j r ) » h = B i m ] - B i u \ + -« liHiQj » ( 78 ) 

j=l 3 

3 1 3 

^( r ) = r CjhimjT) , ?(r ) = — ^ 9j C jH m 3 r o) , j = 1, 2, 3 . 
j=i j'=i 

Before proceeding on to the end of this Section, we discuss the distribution of the roots m| of 
the cubic equation (76) in the complex plane. 

First, ml is always positive since D 3 > 0, as mentioned above, and following the Descartes' 
rule of signs interchange in the sequence of coefficients for real algebraic equations. The other 
two roots are either positive or complex-conjugate with positive real parts. The last 
case leads in (75) to Bessel functions of complex arguments. This fact can indicate that 
the separation of the two groups of functions Pi,( and ©, F r by the ir/2 phase angle, is 
more elaborate than assumed in (57). Another consequence of the existence of complex- 
conjugated roots rrip which is more important from the physical standpoint, is appearance of 
the imaginary contributions in the dispersion equation. This can lead to the complex value 
of the growth rate s = s + ioo, as its solution, and to the non-steady (oscillatory) evolution of 

the jet, e.g. ((z,t) oc <^(ro)e si x e *( w '+ fcz ) ) where to denotes frequency of oscillations. 

4 The freedom to choose the physical parameters of LC seems to admit a degeneration of cubic equation 
(76), when some of the roots mj can coincide in different ways. By virtue, such coincidence is not important, 
since it could occur only at specific wave vectors k* , which the coefficients D^, D\, Dq are dependent upon. 
By the other hand, this kind of degeneration might be interesting if k* is accidentally close to the cut-off wave 
vector kd, when the breakage of the LC jet develops. 



4.2 Dispersion equation 

In what follows we derive the dispersion equation s = s(kr ), which determines the evolution 
of Rayleigh instability in LC jets. The revised version of the boundary conditions (61) at 
r = r , which utilizes the new stream function x( r )> reads 

s JT = k lT> s 7i/^2^- = B 3 — A 2cX + B b k — , -p = 2sp 1 — + T— , (79) 
or or or or or k or z or 

where B 5 = f3 2 + 71^1/^2- Substituting (75) and (78) into (79), and elimination of the coeffi- 
cients Ci, C 2 , C 3 from the linear equations, leads to a (3 x 3)-determinant equation 



det S i:j = , (80) 



where 



hj = k 2 - ^gj, S 2j = B 3 m 2 + B 5 k 2 - ^ 7i/i 2 ^ 



s 3j = r 



1) 



.wij h(mjr Q ) h(mjr 
and /((y) = dli(y)/dy. Equation (80) is an implicit form of the exact dispersion relation, 
which is highly complex and cannot be solved analytically in the general case. Nevertheless, 
here we can verify the fact, that the cut-off wavelength A d does coincide with A s obtained 
due to Plateau. Indeed, the cut-off regime corresponds to (79) when s = and is satisfied 
for T = 0, i.e. A d = A s . The implications of equation (80) can be extended further: for the 
study of different modes of LC flow, including oscillations, and in order to describe asymptotic 
behaviour of LC jets. This is outside the scope of this paper and will be considered elsewhere. 
In the next Section, we present a case which facilitates decoupling of hydrodynamic and 
orientational modes, and consequently the solution of the Rayleigh instability problem in 
closed form. 



5 Decoupling of hydrodynamic and orientational modes 

In this Section we discuss a case that renders the dispersion equation (80) solvable. Here 
we encounter another problem: the elasticity of LC and anisotropy of its viscous properties 
have the same origin and therefore cannot be managed separately. Nevertheless, we consider 
the case where the dispersion equation (80) can be simplified. The large number of physical 



parameters involved (three viscous moduli, two kinetic coefficients A and 71, orientational •& 
and kinematic viscosities z/j, and dimensionless parameter x) call for such a treatment. 
This applies to LC with rod-like molecules (A ~ 1) and low orientational viscosity $ 

A«i ~ 1 , 112 ^ , < Vi , k 2 < I , (82) 

where the first three relations in (82) apply to known nematic LC fluids (see Tables 1, 2 in 
Appendix). The last inequality in (82) applies to the low-viscosity limit which was considered 
for the kinematic viscosity in ordinary liquids by Rayleigh [1]. 
In this case the characteristic equation (76) reduces as follows 

m e --y + -jBe + ^) m *-^(e + ^)=o, = B = h±h. (83) 

The three roots m 2 of equation (76) read 

2m? 2 = Bk 2 + = ± J (B 2 -4)^ + 2(8 -2)k 2 = + ( , m 2 3 = - . (84) 

A simple analysis of (84) shows that the dimensionless parameter B has a critical value 2 
that separates two different evolution scenaria of the LC jet. If B > 2 then the both roots 
m\ and m 2 . are positive and the capillary instability always appears via trivial bifurcation 
(steady-state instability). This scenario applies to MBBA and PA A liquid crystals where 
Bmbba = 5.92, Bpaa = 7.11 (see Tables 1 and 2 in Appendix). In the opposite case, B < 2, 
one can find the regime where the above roots are complex-conjugates. This gives rise to the 
oscillatory evolution of the jet which appears via Hopf bifurcation (see Section 4.1). 

A significant simplification can be obtained if we assume degeneration of the three viscosi- 
ties at critical value B* = 2. Indeed, when the viscous moduli f3j satisfy the relation 

2 — 2/3 1 + /3 3 = 3/3 2 , (85) 

the three roots m 2 of equation (76) are 

m 2 u = k 2 , ml* = k 2 + = , m\ = ^ . (86) 

Note that (85) cancels the last term in (16). The expressions (86) indicate that the problem 
was decomposed in two parts, or, in other words, the cross-terms in equations (70), (71) 



are dropped. Thus, the first part of the problem is associated with Rayleigh instability as 
described by 



(A 2c - ml) (A 2c - ml) X = 0, 



(87) 



with boundary conditions (BC) that account for the elasticity 



d_ 

dr 



dr 



A 2cX + k 2 ^ = , ^-p = 2s(3 1 ^ + T^ at r = r 

k 



d 2 X pdx 
dr 2 dr 



(88) 



The second part is associated with an orientational instability of the director field n(r, t), 



(A 2c - ml) 9 = , with BC s^- = k 2 ^- at r = r . 
\ / dr or 



(89) 



The solutions of equations (87) and (89) are 



X(r) = — /o(mi»r) + —^-Io(m 2 *r) , 9(r) = —I (m 3 r) . 
m u m 2 * m 3 



(90) 



Hence, using these solutions, the hydrodynamic pressure p(r), stream function ip(r) and sur- 
face displacement c(?"o) are obtained as 

p(r) = -cisplQ(m u r) , i/;(r) = r [cj^m^r) + c 2 h(m 2 *r)] , q(r ) = -^h(m 3 r ) , 

where here the only indeterminate are c\ and c 2 , while c 3 can be expressed as their linear 
combination, 



c s_ _ c h(m u r ) ^ h(m 2 *r ) 
k 2 h(m 3 r ) Ji(m 3 r ) 

provided that s = s(kr ) satisfies the dispersion relation which comes from (88), (90) 



(91) 



s 2 + 



" 1:1 r ,,,, v 2&m 2 * i"i(£r ) ,,. 



Io(kr ) 



k 2 +mll 1 (m 2 *r ) 
a ok u , 2 2/1 , An ii(£x ) ml, - /c 2 



(92) 



If x = and z7T = F2, then equation (92) is known as Weber equation for a viscous isotropic 
liquid [6]. For low viscosity, p 1 ~ (3 2 <C y/po^r^, a Rayleigh type expression is obtained (see 
Figure 2) 



(93) 



where subscript "_" denotes low viscosity. 



0.2 0.4 0.6 0.8 1 

k r 



Figure 2: A plot of rescaled growth rate S vs. kr for low viscosity V pr 'q/ctq s_(Ax ) (plain 
line) and high viscosity 2/5 2 r /cr s + (kr ) (dashed line) for different values of x in descending 
order from above: x = 0, 0.25, 1, 5. If $/u — 4x, then the scaling for both viscous regimes 
is the same. 

The maximum s™ ax in equation (93) which corresponds to the wave number k™ ax , gives rise 
to evolution of the largest capillary instability. Numerical calculation shows that s™ aa: and 
fcmax are bokh proportional to (1 + 2x) -1 / 2 



s™ ax S _ 7 == = J-l , fc™ aa; = , a = 0.697 . (94) 



1 pfo_ 
3 VI + 2xV prg ' — r Vl + 2x 

When high viscosity prevails /9i ~ /3 2 ^> y/pd^ro, the dispersion equation reads (see Figure 2) 
a+ (*ro) = ^r ^-^(l + gx)]/^) sr ^^, ^ =0 . (95) 

where subscript " + " denotes high viscosity. Similar to ordinary liquids [4], in this limit there 
is no finite mode of maximum instability for any x. In this case we have 

^) = ^ [cih(k™ ax r ) + c 2 h(m 2 *r )} = . (96) 
s + 

Nevertheless, there exists a continuous range [ 0, (1 + 2x) _1 / 2 r ~ 1 ] of wave numbers k, with 
finite disturbance growth rate s + (kro), which affect the cylindrical jet. 



5.1 Hydrodynamic influence on LC's orientational instability 

We conclude this Section with a brief discussion regarding the hydrodynamic influence on the 
orientational instability of the director field n(r, t). As the effect of hydrodynamics changes 



the wave number k s of Plateau instability to k max , the flow drives the orientational instability 
(43) of the director field n(r, t). Indeed, according to (90) 



nl(r, z) = c 3 h (mr x r) , = ^ — . (97) 

It is convenient to consider the following two marginal viscous regimes. 

1. The low- viscosity limit: 

vTTH? 7S • £ = f • (98) 

where e ~ 1CT 6 -r- 1CT 4 is a small dimensionless parameter. 

2. The high-viscosity limit: 

K+^) 2 -^|- (") 

In both limits the distribution of director field n(r, t) in the jet is always nontrivial and 
definitely far from static distribution (43). 



6 Conclusion 

1. The capillary instability of liquid crystalline (LC) jet is considered in the framework of 
linear hydrodynamics of uniaxial nematic LC. Its static version, called Plateau instability and 
being correspondent to the variational problem of minimal free energy, predicts an essential 
dependence of the disturbance cut-off wavelength upon the dimensionless parameter x = 
K/a r . 

2. The hydrodynamic problem of capillary instability in LC jets is solved exactly followed 
by derivation of the dispersion relation. This relation, which is represented as a determinant 
equation, expresses implicitly the dispersion s = s(k) of the growth rate s, as a function of 
the wave number k of axisymmetric disturbances of the jet. 

3. The case, where the dispersion equation becomes solvable, is considered in detail. 
It corresponds to the regime, wherein the hydrodynamic and orientational modes become 
decoupled. Hydrodynamics changes the wave number k s of Plateau instability into k max that 
produces evolution of the largest capillary instability. Similarly, hydrodynamic flow influences 
the static orientational instability of the director field n(r, t). 



4. The present theory can be easily extended to non-uniaxial nematic LC, which possesses 
finite point symmetry groups G C 0(3) as distinguished from uniaxial group D^. The 
corresponding expressions for the free energy density Ed(G) and the dissipative function D(G) 
were derived in [21]. 

5. In this work the effect of external fields was not considered. However, the theory 
developed here facilitates the treatment of Rayleigh instability in nematic LC in the presence 
of static electromagnetic fields. 
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A Appendix 



Table 1. The basic physical parameters ctj, p, K, a and their derivatives rji, fy, 7$, p h A 
and u,i for nematic liquid crystal ^-metoxybenziliden-^-butilanilin (MBBA) at 25°C taken 

from [19], [20]. 



cui, mPa ■ s 


ai2, 7nPa • s 


o?3, mPa • s 


ai4, mPa • s 


0:5, mPa • s 


ctg, mPa • s 


7 


-78 


-1 


84 


46 


-33 


T)i, mPa • s 


?73, mPa • s 


775, mPa • s 


A 


/ii 




42 


50 


104 


1.026 


1.013 


0.013 


f3i, mPa • s 


#2, mPa • s 


/3 3 , mPa • s 


/?4, mPa • s 


7i, mPa ■ s 


72, mPa • s 


42 


25 


79 


59 


77 


-79 


£>!, mPa ■ s 


P 2) m-P« • s 


B 3 , mPa ■ s 


P 4 , mPa • s 


B 


t?, m 2 /s 


58 


104 


25 


78 


5.92 


1.2 x 10~ 10 


p, kg/m 3 


AT, 


cr , iV/m 


K/<7o, ra 


Vi, m 2 /s 




1.2 x 10 3 


9 x 10~ 12 


38 x 10~ 3 


2.4 x 10~ 10 


1Q- 5 -r- 10~ 4 


10~ 6 -r- 10- 5 



Table 2. The basic physical parameters ccj, p, K, <jq and their derivatives rji, fy, 7$, Sj, //j, A 
and v, L for nematic liquid crystal para-azoxyanisole (PAA) at 122°C taken from [19], [20]. 



cci, mPa • s 


ai2, mPa • s 


CK3, mPa • s 


o?4, mPa • s 


0:5, mPa • s 


ctg, mPa • s 


4 


-6.9 


-0.2 


6.8 


5 


-2.1 


771, mPa • s 


773, mPa • s 


775, mPa • s 


A 


Pi 


P2 


3.4 


4.5 


13.7 


1.06 


1.03 


0.03 


/3i, mPa • s 


P2, mPa ■ s 


/3s, mPa ■ s 


/?4, mPa • s 


71, mPa ■ s 


72, mPa • s 


3.4 


2.25 


11.45 


4.55 


6.7 


-7.1 


Bi, mPa ■ s 


B 2 , mPa ■ s 


B 3 , mPa ■ s 


P 4 , mPa • s 




7?, m 2 /s 


4.34 


9.36 


2.26 


11.24 


7.11 


1.8 x 10~ 9 


p, kg/m 3 


K, N 


(T , 7V/m 


X/(T , 771 


v h m 2 /s 




1.4 x 10 3 


11.9 x 10~ 12 


40 x 10~ 3 


3 x 10~ 10 


10~ 6 -=- 10~ 5 


10~ 4 -=- 10~ 3 
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